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We show tliat tlie middle-range exchange-correlation hybrid of Henderson, Izmaylov, Scuseria 
and Savin (HISS) performs extremely well for elemental and binary semiconductors with narrow or 
visible spectrum band gaps, as well as some wider gap or more ionic systems used commercially. 
The lattice parameters are superior to those predicted by the screened hybrid functional of Heyd, 
Scuseria and Ernzerhof (HSE), and provide a significant improvement over geometries predicted 
by semilocal functionals. HISS also yields band gaps superior to those produced by functionals 
developed specifically for the solid state. Timings indicate that HISS is more computationally 
efficient than HSE, implying that the high quality lattice constants coupled with improved optical 
band gap predictions render HISS a useful adjunct to HSE in the modeling of geometry-sensitive 
semiconductors. 

PACS numbers: 71.15.m 71.15.Mb, 71.20.b, 71.20.Nr, 71.28-|-d, 71.45.Gm, 78.20.-e 



2 



I. INTRODUCTION 



Reliable theoretical treatment of semiconductors, particularly those in the visible range (20 meV-4 eV) remains 
non-trivial. While semilocal DFT methods are the most tractable approach in terms of computational expense, 
the band gaps obtained are often not that impressive. "'^"^ Moreover, many semiconductors, such as the perovskites,^'^ 
chalcopyrites''"^ or Delafossites^^'^"'^ have solid state properties that are very sensitive to geometry. These shortcomings 
of solid state DFT arise, in part, from the various approximations of the exchange-correlation energy, E^c-^^ 

In the simplest case, the Local Spin Density Approximation (LSDA)^'^ assumes that the exchange-correlation 
energy density at the point r depends only on the spin densities at r, rendering the exchange interaction local}'^ For 
solids, the LSDA underestimates lattice constants, ^^'^^ and severely underestimates band gaps."'^'' Generalized gradient 
approximations (GGAs) such as the functional of Perdew, Burke, and Ernzerhof (PBE),^®'^^ improve upon the LSDA 
by incorporating the local density gradient, while the meta-GGAs such as the functional of Tao, Perdew, Staroverov, 
and Scuseria (TPSS)^° also add dependence on the kinetic energy density and/or the Laplacian of the density. While 
GGAs and meta-GGAs are semilocal in character, they can partially account for inhomgeneities in the density, thus 
showing significant improvement over the LSDA in atoms and molecules, but for extended systems, these functionals 
can often provide predictions of lower quality than the LSDA.^^"^'^ 

Hybrid functionals including a percentage of exact non-local Hartree-Fock-type exchange^^'^^ have proven very 
successful for molecules, as is evident from the popularity of the B3LYP^^ functional. These hybrids describe the 
exchange-correlation energy as 

E^^^Eg^^ + c{E^^^E^^^) (1) 

where E^f^ is the exchange-correlation energy of a (semi)local density functional approximation (DFA), E^^ and 
E^^^ are, respectively, the exact exchange energy and the exchange energy from a DFA and, the constant c typically 
ranges from 1/10 to 1/2. Unfortunately, even increasing ly sophisticated hybrids such as PBEh,27^28 TPSSh,^^ or the 
highly empirical MOG"^" series of functionals have been unable to supplant the LSDA for applications in extended 
systems, despite its numerous and well-known insufficiencies. 

In part, this failure is a matter of computational affordability, as the lattice sums required to evaluate the non-local 
exchange energy are slowly convergent, particularly as the band gap of the system decreases. A second problem 
arises from the fact that hybrid functionals are almost always constructed in the generalized Kohn-Sham (GKS) 
framework, '^^ where the the nonlocal exchange energy leads to a non-local exchange "potential" similar to that of 
Hartree-Fock (HF). Accordingly, hybrid functionals formulated in the GKS sense inherit many of the deficiencies that 
Hartree-Fock exhibits for extended systems. In fact, for some solid-state properties, hybrid functionals yield errors 
similar in magnitude to those of semilocal functionals, but at much greater computational cost.'^^''^'^ 

In order to efficiently recover the benefits of non-local exchange for both solids and molecules, a more flexible 
scheme to mix exact and semilocal exchange is required. Probably the most promising route is the use of range- 
separation, ^^'^^"'^^ where the amount of non-local exchange included in the functional depends on the distance between 
electrons. A typical range-separation scheme might split the Coulomb operator l/r into short range (SR) and long 
range (LR) pieces, as 

1 1 — eriiiur) erf(a;r) , , 

- = ^ + ^^, (2) 

r " 




where w is a parameter defining the length scale of separation and erf is the error function, chosen for simplicity of 
integration with gaussian orbitals. The exchange-correlation energy then becomes 

Exc - E^^ + csB [Ex - E^ ) 

so that the amount of non-local exchange used for electrons which are far apart can differ from the amount used for 
electrons which are close together. In molecules, we generally set clr = 1 in what is known as a long-range-corrected 
hybrid. By choosing clr = 1, the correct long-range portion of the exchange-correlation potential is obtained, 
resulting in substantial improvements over non-local hybrids for properties which sample density tails, including non- 
linear optical properties,^'' long-chain polarizabilities,^^ and charge transfer or Rydberg excitations.'^^ Solids, however, 
have no density tails, and the long-range portion of exact exchange becomes problematic if left unscreened, not only 
because it is expensive to compute, but because it yields unphysical results. We, therefore, set c^^ = in what are 
known as screened hybrids. 
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The screened hybrid functional of Hcyd, Scuseria, and Enzcrhof^*^"*^ uses 25% HF exchange in the short range, with 
w = O.lla^"'^, where oq is the Bohr radius. This value of u corresponds to a range of ~ 5 A, so that non-local exchange 
extends over two to three nearest neighbors. Note that for practical purposes this range is quite long. The HSE 
screened hybrid works reasonably well for metals, but significantly improves upon the predictions of regular hybrids 
for semiconductors. For molecular systems, it yields results of the same quality as unscreened hybrids. Nevertheless, 
despite some ability to describe long range properties,"*^ HSE docs not adequately describe quantities sensitive to the 
long-range portion of the exchange-correlation potential. 

A functional successfully coupling the accuracy for solids of a screened hybrid with the accuracy for molecules of 
a long-range-corrected hybrid would be of considerable utility in modeling extended systems, as surface phenomena 
such as molecular adsorption, catalysis, and surface reconstruction could then be reliably addressed. This was the 
motivation behind the development of the HISS functional of Henderson, Izmaylov, Scuseria, and Savin^^^^^ which 
has three ranges, rather than just two. Like HSE, HISS uses FEE as the underlying GGA, but HISS writes 

E^I'' = + CMR {E'J'-"' «-^^^) , (4) 

where the last two terms in parentheses are the middle-range (MR) exact exchange and middle-range FEE exchange 
energies, given by 

1 erfc(wsflr) erf(a;Li?r) evi{LosRr) - evi{ujLRr) 

- = : + 1 + : • (5) 




The three parameters were determined by fitting to small test sets of atomization energies, barrier heights, and 
band gaps, resulting in ujsr. = 0.84aQ *, ujlr. = 0.20aQ *, and cmr = 0.60. The fraction of exact exchange included in 
HISS rises from zero at r = to roughly 0.35 near r = 1.5ao, then decays to zero again for large distances; HISS has 
more exact exchange than does HSE over one or two chemical bonds, but the fraction of exact exchange decays even 
more rapidly than it does in HSE, producing a speed-up in solid-state calculations. 

As expected, enhanced mid-range non-local exchange was demonstrated to reduce calculated errors for thermo- 
chemistry, low-lying Rydberg states, atomization energies, barrier heights and polarizabilities of H2 chains. For the 
semiconductor band gaps tested, HISS did not performs as well as HSE, overestimating them by 0.40 cV.^^ However, 
given the excellent reproduction of molecular properties and the fact that only a small set of solids was examined, we 
chose to more thoroughly assess the performance of HISS in semiconductors. 



II. COMPUTATIONAL METHODS 



Density functional calculations were performed using the periodic boundary-condition code within the GAUS- 
SIAN suite of programs. Gaussian basis sets optimized for solids were discussed and provided in the Supplementary 
Material^* of the earlier study by Heyd and co-workers. 

Eand gaps and related properties reported here are for fully-relaxed, equilibrium structures (lattice parameters 
and geometries) obtained using the LSDA (with SVWNS^^^), FEE, TFSS, FEEsol,^^ the non-local screened hybrid 
functional, HSE, and the middle-range screened hybrid functional HISS. The equilibrium structures are optimized in 
redundant internal coordinates^^ using the GAUSSIAN09 defaults. The SCF convergence criterion was set to 10~^ on 
the density matrix (scf=tight) with Erillouin space integration performed using 12 points in each dimension, for a 
total of 868 points. All extended systems studied are for 3-dimensional "bulk" crystalline materials. 



III. BENCHMARKS 



A. The SC40 Test Set 



The solids used in this study were previously published as the SC40 test sct.^^ The 40 compounds are here separated 
into three groups as presented in Table I. Group 1 contains twenty commercially-relevant bulk semiconductors 
absorbing in the optical (0.20 meV to 4 eV) range. This set includes the most commonly used or benchmarked 
elemental and binary compounds adopting the diamond (di), zincblende (zb), or wurtzite (wu) lattices, and each 
compound is referenced to a complete set of low (LT) and room temperature (RT) band gap data. Group 2 consists 
of a smaller set of wide-gap semiconductors and ionic solids known to be challenging computationally or for various 
intriguing experimental properties. Group 3 encompasses a collection of binary compounds useful for simultaneously 
benchmarking lattice parameters and providing an opportunity for band gap prediction as these materials have gaps 
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that are either unknown or not particularly well-established. Each group is further subdivided according to the lattice 
structure, with the diamond, zincblende and wurtzite structures preceding any rock salt (rs) structures within each 
group. 

To facilitate evaluation of the relative performance of each functional, the experimental band gaps have been updated 
to reflect half a decade of experimental progress and further augmented to include, and discriminate between, room 
temperature (RT) and low temperature (LT) values. It should be noted that the low temperature data set is not 
exhaustive and the experimental measurements range from extrapolated K values up through 27 K, providing a 
rather inconsistent basis by which to judge the quality of the predicted optical gaps, particularly if the temperature 
dependence of a given material is uncertain. Alternatively, the more complete room temperature (RT) data set, 
which has gaps typically ca. 0.2 eV higher in energy than the measured LT results, can provide a useful baseline 
for estimations of the low temperature gaps where unavailable or extrapolated. The majority of the experimental 
reference data (gaps and lattice constants) are taken from the 3rd Edition of the Springer Semiconductor Handbook^^ 
and any references included therein. Alternate tabulations are cited where applicable. 



B. Lattice Parameters 

The LSDA provides reasonable, but generally contracted lattices, while gradient-corrected functional tend to 
over-correct and expand them, instigating the development of functionals specifically-tailored to produce reasonable 
solid state geometries. Arguably, the most well-known of these special functionals is PBEsol,'"'^ a modified GGA 
with a significantly improved lattice constants for solids (^4x better) relative to PBE. Nevertheless, the virtues 
of PBEsol were not seen to extend to semiconductors.^^ Interestingly, a recent assessment of meta-GGAs,^'' for 
twenty solids utilizing experimental lattice parameters extrapolated to K and then corrected with the Zero Point 
Anharmonicity Expansion (ZPAE), indicates that LSDA does not perform as well as any of the special functionals 
studied. Nevertheless, markedly different trends in performance arc seen upon examination of a subset including 
only the five semiconductors and MgO present in the aforementioned test set. A comparison of these errors, relative 
to corrected low temperature experiment, for LSDA, PBEsol, AM05,^® and revTPPS,^^ is extended to include the 
screened exchange hybrids HSE and HISS, as is summarized in Table II. 

For this test set, the mean error (ME) of LSDA is actually larger than that of PBEsol for the corrected low 
temperature lattice parameters, and of similar magnitude (but opposite sign) relative to revTPSS. HSE has the 
largest ME, but HISS - not LSDA - provides the best lattice parameters. Considering the mean absolute error 
(MAE) instead, LSDA is only 0.1 pm better than PBEsol, while HISS again proves superior to both. Note also that 
the middle-range screened hybrid yields better lattice parameters for the narrow- and wide-gap cubic semiconductors, 
as well the ionic MgO, even for the uncorrected extrapolations to absolute zero. 

Admittedly, this test set of six solids is rather small, and while the earlier semiconductor test set utilized for HISS^^ 
showed promising results for a slightly different set from that contained in Table II, it is also small. We thus examined 
the performance of HISS for the entire SC40 test set and compared the data to that produced by the LSDA, TPSS, 
PBEsol, TPSS, and HSE functionals. This provides a much larger semiconductor test set while allowing for analysis 
of some wider-gap materials commonly used in conjunction with the most commonly used materials in Group 1 of 
Table I. Given that the PBEsol lattice constants are of higher quality than PBE, we do not include data for the latter 
in this table for aesthetic reasons. Nevertheless, since both GGAs yield band gaps of similar magnitude, PBE data is 
presented in other analyses, particularly those in Section IIIC. 

To better illustrate the magnitudes and systematics of lattice parameter predictions. Figure 1 presents the errors for 
all forty-three lattice parameters contained in Table I, as well as those for PBE, in the form of histograms. The ME for 
each functional is depicted by the vertical dotted line on each plot. The results from HISS are sharply peaked around 
zero error, it is clear that the LSDA underestimates lattice constants, the well-known PBE and TPSS overestimations 
are evident, and HSE produces highly-clustered overestimations with accuracy comparable to that of the LSDA. The 
PBEsol distribution is somewhat narrower than that of LSDA, but, as with other GGAs, is slightly more skewed 
toward expanded lattice constants, with the notable exception of the Ca, Sr and Ba chalcogenides. PBEsol actually 
contracts these Group 3 lattices and performs more like LSDA, particularly in the case of BaTe, the obvious outlier 
in Figure 1. 

As was demonstrated from the smaller test set in Table II, it is helpful to consider both ME and MAE in performance 
evaluations. The mean absolute errors of the equilibrium lattice parameters for the twenty commercially-relevant 
semiconductors in Group 1 are thus presented in Figure 2. 

For this particular group, PBEsol does outperform LSDA, but not HISS. Remarkably, HISS provides excellent 
lattice parameters, with a MAE of less than half that of LSDA, (see Table III) despite producing bond lengths in 
molecules that are not particularly impressive. In fact, when comparing the ten largest deviations from experiment, 
HISS outperforms all semilocal functionals other than PBEsol for the zinc and cadmium chalcogenides, produces more 
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TABLE I. Equilibrium lattice parameters (A) and band gaps (eV) at the relaxed geometries for the SC40 solid test set.^^'^'^ 
Initial lattice constants are room temperature experimental values except for the six presented in Table II. No ZPAE corrections 
have been made. The reference set for band gaps is either low temperature (LT) or room temperature (RT). 



Expt. Lattice'' 



Calculated Lattice Parameters 



Expt. Gap"" 



Calculated Band Gaps 



Solid 


Type ao or a, c 


HISS 


HSE 


LSDA 


TPSS 


PBEsol 


LT 


300 K 


Type 


HISS 


HSE 


LSDA 


TPSS 


PBEsol 


Group I: Semiconductors 


























Ge 


di 


5.658 


5.660 


5.706 


5.650 


5.743 


5.693 


0.74 


0.66 




1.08 


0.54 


0.00 


0.00 


0.00 


Si 


di 


5.430 


5.430 


5.451 


5.410 


5.465 


5.442 


1.17 


1.12 




1.45 


1.22 


0.51 


0.74 


0.53 


;3-SiC 


di 


4.358 


4.352 


4.376 


4.355 


4.395 


4.380 


2.42 


2.20 




2.74 


2.32 


1.33 


1.35 


1.27 


BP 


zb 


4.538 


4.527 


4.547 


4.510 


4.566 


4.540 


2.40 


2.20 




2.43 


2.12 


1.28 


1.41 


1.24 


AlP 


zb 


5.463 


5.459 


5.481 


5.438 


5.498 


5.472 


2.50 


2.45 




2.71 


2.44 


1.53 


1.83 


1.56 


AlAs 


zb 


5.661 


5.668 


5.697 


5.639 


5.713 


5.682 


2.23 


2.15 




2.40 


2.16 


1.33 


1.64 


1.37 


GaP 


zb 


5.451 


5.457 


5.490 


5.419 


5.523 


5.468 


2.35 


2.27 




2.67 


2.42 


1.53 


1.90 


1.62 


GaAs 


zb 


5.648 


5.668 


5.711 


5.626 


5.745 


5.687 


1.52 


1.42 


D 


1.86 


1.18 


0.43 


0.52 


0.42 


GaSb 


zb 


6.096 


6.099 


6.146 


6.043 


6.182 


6.111 


0.82 


0.73 


D 


1.31 


0.70 


0.09 


0.09 


0.06 


InP 


zb 


5.869 


5.872 


5.914 


5.840 


5.961 


5.891 


1.42 


1.34 


D 


2.23 


1.61 


0.83 


0.91 


0.83 


InAs 


zb 


6.058 


6.075 


6.125 


6.038 


6.172 


6.099 


0.42 


0.35 


D 


0.93 


0.36 


0.00 


0.00 


0.00 


InSb 


zb 


6.479 


6.486 


6.04U 


c A on 


0.080 


6.0U1 


U.24 


n 1 o 
O.lo 


U 


U.oO 


U.2o 


U.OO 


U.OO 


n nn 
O.UU 


ZnS 


zb 


5.409 


5.413 


0.40D 


o.oiy 


^ AdK 

0.4D0 


O.ooo 


Q fin 


3.66 


U 


A TO 
4. iZ 


Q Q7 


Z.ZD 


AC\ 
Z.4U 


z.zz 


ZnSe 


zb 


5.668 


5.686 


5.713 


5.589 


5.738 


5.657 


2.82 


2.50 


U 


2.93 


2.27 


1.21 


1.48 


1.26 


ZnTe 


zb 


6.089 


6.127 


a 1 

6.iOD 


6.U1 1 


0.1 (0 


G non 

6.uoy 


o on 
2.39 


2.30 


L) 


Z.I / 


z.ib 


1 oo 


i.4o 


1 on 


CdSe 


zb 


6.052 


6.118 


6.156 


6.024 


6.195 


6.098 


1.90 


1.74 


U 


1.90 


1.36 


0.34 


0.73 


0.45 


CdTe 


zb 


6.480 


6.534 


6.575 


6.423 


6.611 


6.502 


1.57 


1.48 


U 


2.00 


1.49 


0.61 


0.88 


0.67 


/?-GaN 


wu 


4.523 


4.483 


4.520 


4.476 


4.552 


4.520 


3.30 


3.23 


U 


4.05 


2.97 


1.93 


1.56 


1.70 


a-GaN 


wu (a) 


3.189 


3.176 


3.197 


3.164 


3.222 


3.197 


3.50 


3.44 


U 


4.23 


3.14 


2.08 


1.74 


1.85 




wu (b) 


5.185 


5.167 


5.215 


5.171 


5.246 


5.216 


















InN 


wu (a) 


3.537 


3.524 


3.555 


3.518 


3.585 


3.552 


0.72 


0.69'' 


D 


1.51 


0.66 


0.00 


0.00 


0.00 




wu (c) 


5.704 


5.682 


0. 1 4U 


o.oy ( 


0. / / D 


0. ( 40 


















Group II: Wide 


-Gap/Notable 












5.48 














C 


di 


3.567 


3.536 


3.556 


3.538 


3.580 


3.563 


5.50 


I 


6.11 


5.42 


4.19 


4.19 


4.05 


AlSb 


zb 


6.136 


6.138 


6.172 


6.098 


6.194 


6.146 


1.69 


1.62 


I 


2.05 


1.85 


1.18 


1.49 


1.22 


BN 


zb 


3.616 


3.586 


3.605 


3.585 


3.629 


3.611 


6.36 


6.20 


I 


6.69 


5.90 


4.40 


4.47 


4.30 


CdS 


zb 


5.818 


5.865 


5.899 


5.776 


5.944 


5.844 


2.55 


2.31 


D 


2.72 


2.10 


1.01 


1.34 


1.08 


MgS 


zb 


5.622 


5.670 


5.690 


5.619 


5.518 


5.675 


4.78 




D 


5.17 


4.48 


3.37 


3.68 


3.34 


AIN 


wu (a) 


3.111 


3.109 


3.124 


3.105 


3.141 


3.131 


3.60 


3.50 


I 


6.62 


5.50 


4.29 


3.96 


3.97 




wu (c) 


4.981 


4.967 


4.999 


4.968 


5.030 


5.009 


















MgTe 


zb 


6.420 


6.445 


6.469 


6.362 


6.504 


6.432 


6.19 


6.13 


D 


3.91 


3.49 


2.57 


2.96 


2.58 


MgO 


rs 


4.207 


4.200 


4.222 


4.178 


4.249 


4.231 


7.90 


7.30 


D 


7.87 


6.40 


4.90 


4.58 


4.50 


Group III: Benchmarking'^ 


























BAs 


zb 


4.777 


4.769 


4.797 


4.750 


4.820 


4.788 


— 


1.46 


I 


2.14 


1.89 


1.11 


1.24 


1.10 


MgSe 


rs 


5.460 


5.484 


5.508 


5.418 


5.518 


5.476 




2.47 




3.05 


2.58 


1.62 


2.01 


1.69 


CaS 


rs 


5.689 


5.696 


5.704 


5.582 


5.710 


5.640 








4.02 


3.52 


2.10 


2.62 


2.27 


CaSe 


rs 


5.916 


5.946 


5.954 


5.801 


5.954 


5.879 








3.39 


2.97 


1.68 


2.17 


1.85 


CaTe 


rs 


6.348 


6.385 


6.390 


6.211 


6.388 


6.294 








2.66 


2.32 


1.19 


1.67 


1.35 


SrS 


rs 


5.990 


6.030 


6.043 


5.924 


6.055 


5.985 




4.10 




3.97 


3.52 


2.20 


2.73 


2.38 


SrSe 


rs 


6.234 


6.277 


6.289 


6.159 


6.290 


6.218 








3.41 


3.03 


1.85 


2.31 


1.99 


SrTe 


rs 


6.640 


6.698 


6.707 


6.539 


6.706 


6.615 








2.82 


2.51 


1.42 


1.93 


1.59 


BaS 


rs 


6.389 


6.401 


6.414 


6.298 


6.433 


6.350 




3.88 




3.61 


3.21 


1.99 


2.55 


2.15 


BaSe 


rs 


6.595 


6.640 


6.653 


6.515 


6.658 


6.577 




3.58 




3.14 


2.80 


1.69 


2.18 


1.83 


BaTe 


rs 


7.007 


7.002 


7.018 


6.854 


7.024 


6.919 




3.08 




2.48 


2.22 


1.25 


1.73 


1.38 


BSb 


zb 


5.278 


5.222 


5.256 


5.201 


5.280 


5.240 






D 


1.48 


1.31 


0.74 


0.76 


0.70 



^ The abbreviations di, zb, wu and rs refer to the diamond, zincblende, wurtzite and rock salt lattices, respectively. The wurtzite 
structure is hexagonal, so there are two lattice vectors a and c. 

^ The abbreviations I or D refer to indirect or direct gap transitions; all experimental data taken from Reference 56 and references 
therein, unless otherwise noted. Low temperature experimental band gaps vary in terms of temperature, and are not available for all 
systems: K (extrapolated): Si, GaP, GaAs, GaSb, CdTe^S, MgS^O, MgO; < 2K: Ge, /3-SiC, AlP, InP, InSb, ZnTe, a-GaN, 
MgTe^l; ~ 4 - 6K: AlAs, InN,'^^ j^^g^ gnSc; ~ 7 - lOK: /3-GaN, BN,''^ AIN; 14 K: ZnS''-* ; 27 K: AlSb ; Others: CdSe,'^^ CdS,''^ 
BAs. 

All experimental data from Table I, Reference 52 except for MgSc (Reference 56.) 
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accurate lattice parameters for the a-GaN and InN wurtzite structures and excels for the diamond lattices. HISS, 
like PBEsol also performs well for Groups 2 and 3, while the LSDA errors are notably larger for Group 3. All of the 
error data is summarized by group in Table III. 



TABLE II. Comparison of equilibrium lattice parameters for five semiconductors and MgO relative to uncorrected and 
corrected extrapolations to K from experimental values, as provided in Ref. 67. Mean error (ME) and mean average errors 
(MAE) reported in pm. (Corrected data is obtained via subtraction of the zero-point anharmonic expansion (ZPAE) from the 
extrapolated K values. HISS and HSE values are from this work, the remainder are from Ref 67.) 









Theoretical 


ao (A) 






Experimental ao (A) 


Solid 


LSDA 


PBEsol 


AM05 


revTPSS 


HSE 


HISS 


K-ZPAE 


K 


Ge 


5.631 


5.680 


5.685 


5.682 


5.706 


5.660 


5.644 


5.658 


Si 


5.405 


5.432 


5.434 


5.439 


5.451 


5.430 


5.421 


5.430 


P-SiC 


4.332 


4.355 


4.352 


4.357 


4.376 


4.352 


4.346 


4.358 


GaAs 


5.615 


5.665 


5.672 


5.680 


5.711 


5.668 


5.640 


5.648 


C 


3.533 


3.552 


3.548 


3.558 


3.556 


3.535 


3.553 


3.567 


MgO 
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4.222 


4.200 
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4.207 


ME 


-2.5 


1.9 
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3.8 


0.9 


ZPAE-corrected 
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1.9 


2.3 
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3.8 
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extrapolation to K 




ME 


-3.0 


0.7 
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1.5 


2.6 


-0.4 


Uncorrected 




MAE 


3.0 


1.3 


1.7 


1.8 


2.9 


1.1 


extrapolation to K 
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TPSS 



PBEsol 



i 
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XL 
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FIG. 1. Histograms of calculated lattice constants errors (pm) relative to experimental values from Table I. The dashed lines 
indicate MEs for HISS, HSE, LSDA, PBE, TPSS and PBEsol: 0.9, 3.7, -4.8, 7.3, 6.0 and 0.5 pm, respectively. 
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TABLE III. Error Statistics for equilbrium geometries (pm) and band gaps (eV) calculated at those geometries. The band gap 
data is compared relative to low (LT) and room (RT) temperature experiment. Groups 1-3 defined in Table I. 



Lattice Parameters (pm) Band Gaps (eV) 
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3.41 
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2.09 


MAE 


0.39 


0.21 


0.98 


0.88 


1.00 


Groups 1+2 














LT Groups 14-2 














ME 


0.50 


3.74 


-3.22 


6.71 


1.75 


ME 


0.38 


-0.24 


-1.15 


-0.93 


-1.19 




MAE 


1.87 


3.90 


3.22 


6.71 


2.07 


MAE 


0.38 


0.28 


1.15 


0.93 


1.19 


Group 3 














RT Group 1 














ME 


1.88 


3.42 


-8.93 


4.27 


-2.85 


ME 


0.50 


-0.07 


-0.88 


-0.78 


-0.89 




MAE 


3.05 


3.79 


8.93 


4.27 


3.30 


MAE 


0.50 


0.15 


0.88 


0.78 


0.89 


All Groups 














RT Groups 1+2 














ME 


0.89 


3.65 


-4.81 


6.03 


0.46 


ME 


0.49 


-0.12 


-1.02 


-0.93 


-1.06 




MAE 


2.20 


3.87 


4.81 


6.03 


2.42 


MAE 


0.49 


0.20 


1.02 


0.93 


1.06 



^ MAE=Mean Absolute Error; Max=Maximum absolute error 




HISS PBEsol LSDA HSE TPSS PEE 



FIG. 2. MAEs comparing equilibrium/relaxed lattice parameters to the experimental data listed in Group 1 of Table I. This 
data is summarized in Table III. Notice the exceptional performance of HISS. 



C. Band Gaps 



It should be emphasized that screened hybrids, such as HISS and HSE, which calculate nonlocal Hartree-Fock 
type exchange with the GKS approach, predict fundamental gaps as band energy (conduction band - valence band) 
differences. ^° They also provide excellent approximations of optical gaps, particularly for visible-range semiconductors 
and other materials. ^-"^ This is well-illustrated in Figure 3, which plots the mean absolute errors (MAEs) for the 
calculated HISS, HSE, LSDA and PBEsol band gaps of Group 1 (Tabic I) against low temperature experimental gaps. 
The dashed lines indicate deviations of ±0.3 cV. 

Almost all HSE data is within the ±0.3 eV range, while HISS tends to run toward larger gaps and the semilocal 
functionals consistently underestimate all band gaps by more than 0.3 eV. Our group previously demonstrated that 
HSE is capable of correctly describing the Mott metal-insulator transition in actinide oxides: materials known to be 
narrow-gap semiconductors. ^^"^^ All of the semilocal functionals predict these systems to be metallic. Similarly, for 
the SC40 test set, the LSDA metallizes the narrow band gaps of Gc and the indium pnictides as is evident from he 
clustering on the x-axis between and 1 in Figure 3. HSE - and HISS - not only predict narrow gaps, but do so in 
good agreement with experiment. 

It is interesting to note that while PBEsol and HSE yield similar lattice parameters for nitride compounds adopting 
the wurtzite structure, the band gaps are significantly different: HSE band gaps are in much better agreement 
with experiment. As is evident from Table I, HISS consistently produces better lattice parameters than semilocal 
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Experimental Gap (eV) 

FIG. 3. Comparison of DFT band gaps vs. low temperature experimental gaps, for Group I of Table I. The 20 semiconductors 
have gaps in the visible range. Dotted lines indicate deviations of ± 0.3 eV. 



functionals accompanied by systematically-larger band gaps. These overestimations are, nevertheless, typically smaller 
in magnitude than the well-known undcrestimations expected from semilocal functionals. 

The MAEs for Groups 1 and 2 are summarized on the right of Table III. As expected, HSE shows the smallest 
errors, on the order of 0.20 eV for either temperature data set, but HISS deviates less for more ionic compounds like 
MgO or the Mg chalcogenides, and yields better lattice parameters. The performance of all functionals for Group 1 
semiconductors relative to both temperature ranges is summarized in in Figure 4. In marked contrast to all other 
functionals, HISS appears to have errors in better agreement with respect to the reported LT gaps than the RT gaps, 
an artifact of the systematic overestimation of band gaps demonstrated in Figure 3. 




LSDA PBE PBEsol TPSS HSE HISS 



FIG. 4. Mean absolute errors relative to both low (gray) and room (blue) temperature experimental band gaps for the Group 
Isemiconductors in Table I. The 20 semiconductors absorb in the visible range and adopt diamond, zincblende and wurtzite 
structures. 
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D. Noteworthy Cases 



Having discussed the bulk of the statistics above, we now turn our attention to a few systems worthy of comment. 
As is pointed out in the Springer Handbook, there is conflicting experimental data regarding the magnitude and, at 
one time, some uncertainty regarding the location of the AlSb indirect transition. One set of experiments indicates 
that the indirect transition occurs somewhere near F — > X, with a magnitude of 1.6-1.7 eV,'^^'''^ while another set 
supports a r — !• L transition^'' with a gap closer to 2 cV, and in recent texts, the gap is reported to be ca. 1.5, 1.6 or 
1.7 eV in Refs. 78 and 79 and 80, respectively. 

We calculated the band structure of AlSb at the equilibrium geometries produced by LSDA, HSE and HISS. All 
functionals indicate that the indirect transition is clearly in the vicinity of near F — X, Figure 5. The K HSE gap of 
1.82 eV compares favorably with the 1.69 eV (27 K) experimental band gap, while showing the characteristic camel's 
back structure. 

The slightly larger 2.05 eV band gap predicted by HISS might encourage preference of the F L transition, but as 
is seen from plot to the right of Figure 5, the smallest indirect gap again occurs near X and there is little difference 
between the structure of highest valence and lowest conduction bands for either screened hybrid calculation. The fact 
that HISS systematically overestimates the band gaps, yet shows similar structure to HSE also tends to favor the 
smaller transition at X. Note also that the width of the valence band predicted by HISS is intermediate between that 
of LSDA and HSE. 




U W K I U X 



U W K I U X 



FIG. 5. Band structures for the relaxed geometries of AlSb as calculated with LSDA (left), HSE (middle) and HISS (right). 
Note that all functionals show that the two observed experimental transitions near T ^ X and F — i- L are close in energy, but 
that the F — >■ X gap is smaller. 



The zincblcndc polymorph of CdS, Hawleyite, included in the SC40 test set is not as widely used as wurtzite CdS, 
Greenockite: the hcxagaonal form is commonly employed in thin film solar cells and quantum dots. Note that most 
of the available, reliable experimental data for the cubic CdS is for thin films, which typically have larger band gaps 
than for the bulk, which is what we calculate here. Nevertheless, the HISS and HSE results are within a reasonable 
range and the consistently lower errors HISS yields compared to those of scmilocal functionals are encouraging. 

Finally, the wide-gap nitrides of aluminum and boron are worthy of comment because BN is an unsual zincblende 
compound, while the isoelectronic AIN adopts a rather polar wurtzite structure. Boron nitride, "borazone" is almost 
as hard as diamond, has excellent refractory properties,^''' and, atypical of most zincblende structures, has a large band 
gap (6.3 eV), well into the UV spectrum. Interestingly, the predicted screened hybrid band gaps for these systems 
(Table I) remain in good agreement with experiment, implying some degree of utility for higher-energy semiconductor 
applications. 



E. Timings 



In addition to excellent performance for lattice parameters, HISS possesses an another advantage over HSE: because 
the range-separation parameter lulr for HISS is roughly twice the parameter a; of HSE, the long-range exchange in 
HISS is screened more aggressively. Accordingly, the computational cost for HISS is smaller than that for HSE. The 
claim in Reference 81 that HISS is quite expensive because one must do the two-electron integrals twice is simply 
incorrect. While this is true in the most naive implementation possible, computationally efficient implementations 
of HISS will calculate only the so-called [0]*^™^ integrals twice, and the expensive steps of contraction and angular 
moment incrementation are performed once, just as in HSE. 

Table IV contains the CPU time per SCF cycle relative to that of PBE for BP, C, InN, MgO, Si, and SiC. This set 
includes both wide gap and narrow gap materials adopting diamond, zincblende or wurtzite lattices. All calculations 
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arc performed at the experimental geometries. It is clear that HISS is only twice as expensive as PBE, and that HSE 



TABLE IV. CPU time per SCF cycle relative to PBE for representative systems. The numbers in parentheses are the number 
of SCF cycles required for convergence. 



Solid 




Relative Avg, 


. Time (# of SCF Cycles) 




PBE 


HISS 


HSE06 


PBEh 


C 


(9) 


1.95 (9) 


2.63 (10) 




Si 


(10) 


1.70 (12) 


2.38 (10) 


33.18 (11) 


BP 


(8) 


1.88 (9) 


2.56 (9) 


37.75 (9) 


SiC 


(8) 


1.89 (8) 


2.64 (9) 


39.51 (9) 


MgO 


(10) 


1.85 (10) 


2.73 (8) 


32.77 (15) 


InN 


(15) 


2.21 (9) 


2.91 (7) 




Mean 


1 


1.91 


2.64 


35.50 



costs about 50% more than HISS, both producing much better band gaps and lattice parameters. We also include 
results for the PBEh^^ global hybrid (also known by the acronyms PBEO and PBEIPBE) to illustrate the enormous 
savings resulting from range separation. PBEh calculations on C and InN were not possible, nevertheless, the PBEh 
band gaps for Si, SiC BP, and MgO are respectively, 1.81, 3.22, 2.74, and 7.31 eV. The PBEh band gaps are in error 
by approximately 0.7 eV, offering minimal improvement over PBE but at an enormously larger (40 x) computational 
cost. 

Global hybrid functionals like B3LYP and PBEh have been benchmarked against HSE for a few systems and 
properties, and it has been shown that their increased computational cost (observed for both with plane waves and 
gaussian orbitals) does not lead to improved answers compared to HSE. This is even more true for semiconductors: 
both B3LYP and PBEh overestimate the band gaps*^ 



IV. CONCLUDING REMARKS 

The middle-range screening utilized in the HISS functional allows for accurate reproduction of lattice parameters 
with concomitant band gap prediction at nearly the quality of HSE. The HISS screened hybrid band gap errors 
are systematically over-estimated, but are of increased accuracy for larger gap systems, while lattice parameters 
arc superior to those produced by LSDA - or even PBEsol - for commercially-relevant semiconductors, including 
non-cubic systems. 

In fact, as was recently demonstrated for the perovskite SrTiOa system, ^'^ HISS can rapidly reproduce difficult 
lattice distortions while producing near-HSE-quality band gaps. The predictive power provided by such systematic 
errors was also nicely demonstrated by the disputation of the experimental geometry for the novel Delafossite material 
CuB02,®^"^® where both the accuracy of the band gaps and the generally systematic deviations in lattice parameters 
produced by HSE were used to justify theoretical arguments. Given that HISS performs well for molecules, it is 
reasonable to assume that molecule-surface interactions will show an improvement over what is currently available. 

Finally, timing data indicates that HISS is more computationally efficient than HSE, yet more accurate than 
semilocal functionals, PBEsol or PBEh, imbuing HISS with yet another advantage: reliable theoretical data for 
semiconductors may be acquired by running a calculation that will complete at a computational expense of merely 
twice that of PBE, without necessitating additional high-level corrections. Such easily-obtained, high quality lattice 
constants coupled with improved optical band gap predictions will render HISS a useful adjunct to HSE in theo- 
retical investigations of geometry-sensitive semiconductors and many wider-gap materials requiring accurate lattice 
parameters in order to obtain reliable solid state properties. 



ACKNOWLEDGMENTS 

The work at Rice University was supported by the Qatar National Research Fund (QNRF) through the National 
Priorities Research Program (NPRP Grant No. 09-431-076) and The Welch Foundation (C-0036). The authors would 
like to thank Dr. Donghyung Lee and Mr. Ireneusz W. Bulik helpful discussions and technical assistance. 



11 



G. E. Scuseria and V. N. Staroverov, in Theory and Applications of Computational Chemistry: The First 40 Years, edited 
by C. E. Dykstra, G. Frenking, K. S. Kim, and G. E. Scuseria (Elsevier, Amsterdam, The Netherlands, 2005) pp. 669-724. 
A. Dreuw and M. Head-Gordon, Chem. Rev., 105, 4009 (2006). 

Y. Zhao and D. G. Truhlar, J. Phys. Chem. A, 110, 13126 (2006). 

M. J. G. Peach, P. Benfield, T. Helgaker, and D. J. Tozer, J. Chem. Phys., 128, 044118 (2008). 
M. Grundmann, The Physics of Semiconductors (Springer- Verlag, Berlin, 2010) Chap. 14, p. 427. 
N. lies, F. Finocchi, and K. D. Khodja, Journal of Physics: Condensed Matter, 22, 305001 (2010). 
Wide-Gap Chalcopyrites (Springer, Berlin, 2006). 
Ternary Chalcopyrites (Pergamon, Oxford, 1975). 

M. Grundmann, The Physics of Semiconductors (Springer- Verlag, Berlin, 2010) Chap. 3, pp. 51-53. 
R. G. Gordon, MRS BuU., 8, 52 (2000). 

M. Grundmann, The Physics of Semiconductors (Springer- Verlag, Berlin, 2010) Chap. 19, pp. 511-515. 

J. P. Perdew and K. Schmidt, in Density Functional Theory and Its Applications to Materials, edited by V. V. et al. 

(American Institute of Physics, New York, 2001). 

J. P. Perdew and Y. Wang, Phys. Rev. B, 45, 13244 (1992). 

J. P. Perdew and K. Schmidt (American Institute of Physics, 2001). 

A. Khein, D. J. Singh, and C. Umrigar, Phys. Rev. B, 51, 4105 (1995). 

S. Kurth, J. P. Perdew, and P. Blaha, Int. J. Quantum Chem., 75, 889 (1999). 

J. P. Perdew, Int. J. Quantum Chem. Symp, 19, 497 (1986). 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett., 77, 3865 (1996), E78, 1396 (1997), Ref.l9. 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett., 78, 1396 (1997). 

J. Tao, J. P. Perdew, V. N. Staroverov, and G. E. Scuseria, Phys. Rev. Lett., 91, 146401 (2003). 

S.-K. Ma and K. A. Brueckner, Phys. Rev., 165, 18 (1968). 

J. P. Perdew and Y. Wang, Phys. Rev. B, 33, 8800 (1986), 40, 3399(E) (1986). 

J. P. Perdew, Phys. Rev .B, 33, 8822 (1986), 34, 7406(E) (1986). 

A. D. Becke, J. Chem. Phys., 98, 1372 (1993). 

A. D. Becke, J. Chem. Phys., 98, 5648 (1993). 

P. J. Stephens, F. J. Devlin, C. F. Chabalowski, and M. J. Frisch, J. Phys. Chem., 11623 (1994). 

C. Adamo and V. Barone, J. Chem. Phys., 110, 6158 (1999). 

M. Ernzerhof and G. E. Scuseria, J. Chem. Phys., 110, 5029 (1999). 

V. N. Staroverov, G. E. Scuseria, J. Tao, and J. P. Perdew, J. Chem. Phys., 119, 12129 (2003). 
Y. Zhao and D. G. Truhlar, Theor. Chem. Ace, 120, 215 (2008). 

A. Seidl, A. Gorling, P. Vogl, J. A. Majewski, and M. Levy, Phys. Rev. B, 53, 3764 (1996). 
T. M. Henderson, B. G. Janesko, and G. E. Scuseria, J. Phys. Chem. A, 112, 12530 (2008). 
J. Paier, M. Marsman, G. Kresse, I. C. Gerber, and J. G. Angyan, J. Chem. Phys., 124, 154709 (2006). 
A. Savin, in Recent Developments and Applications of Modern Density Functional Theory, edited by J. M. Seminario (Else- 
vier, Amsterdam, The Netherlands, 1996) pp. 327-357. 
A. Savin and H.-J. Flad, Int. J. Quantum Chem., 56, 327 (1995). 

T. Leininger, H. StoU, H.-J. Werner, and A. Savin, Chem. Phys. Lett., 275, 151 (1997). 
O. A. Vydrov and G. E. Scuseria, J. Chem. Phys., 125, 234109 (2006). 

H. likura, T. Tsuneda, T. Yanai, and K. Hirao, J. Chem. Phys., 115, 3540 (2001). 
H. Sekino, Y. Maeda, and M. Kamiya, Mol. Phys., 103, 2183 (2005). 

J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys., 118, 8207 (2003). 

A. V. Krukau, O. A. Vydrov, A. F. Izmaylov, and G. E. Scuseria, J. Chem. Phys., 125, 224106 (2006). 

T. M. Henderson, A. F. Izmaylov, G. Scalmani, and G. E. Scuseria, J. Chem. Phys., 131, 044108 (2009), specifically, the 
HSEh parameterization of HSE06''^'**'', called by the GAUSSIAN keyword HSEhlPBE. 

B. G. Janesko, T. M. Henderson, and G. E. Scuseria, Phys. Chem. Chem. Phys., 11, 443 (2009). 
V. Barone, O. Hod, J. E. Peraha, and G. E. Scuseria, Acc. Chem. Res., 44, 269 (2011). 

T. M. Henderson, A. F. Izmaylov, G. E. Scuseria, and A. Savin, J. Chem. Phys., 127, 221103 (2007). 

T. M. Henderson, A. F. Izmaylov, G. E. Scuseria, and A. Savin, J. Theor. Comput. Chem., 4, 1254 (2008). 

K. N. Kudin and G. E. Scuseria, Phys. Rev. B, 61, 16440 (2000). 

K. N. Kudin and G. E. Scuseria, Chem. Phys. Lett., 289, 611 (1998). 

K. N. Kudin and G. E. Scuseria, Chem. Phys. Lett., 283, 61 (1998). 

M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, 
B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, 
J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, 
H. Nakai, T. Vreven, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, 
V. N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, 



12 



62 
63 



64 
65 



N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, 
O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. 
Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, O. Parkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, and 
D. J. Fox, "Gaussian 09 Revision B.l and gdv-hll," (2009), Gaussian Inc., Wallingford CT. 

The EPAPS Document No. E-JCPSA6-123-301539 contains modified basis sets, input files and geometries relevant to the 
material in Ref. 52. Please see doi : 10 . 1063/1 . 2085170 for the direct link or http://www.aip.org/pubservs/epaps.html 
for more information on EPAPS. 

J. Heyd, J. E. Peralta, G. E. Scuseria, and R. L. Martin, J. Chem. Phys., 123, 174101 (2005). 
S. H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys, 58, 1200 (1980). 



J . P. Perdew, A. Ruzsinszky, G . I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Phys 

R |ev. Lett., 100, 136406 (2008)| , see also the 2009 Erratum in Reference**. 

K. N. Kudin, G. E. Scuseria, and H. B. Schlegel, J. Chem. Phys., 114, 2919 (2001). 

O. Madelung, Semiconductors - Data Handbook 3rd Ed. (Springer, New York, 2004). 

B. Arnaud qv. T. Paskova. P. P. Paskov. B. Magnusson, E. Valcheva, B. Monemar, H. Lu, W. J. SchafT, H. Amano, and 



I. Akasaki, [Phys. Rev. B, 69, 115216 (2004) 



B. G. Yacobi, Semiconductor Materials: An Introduction to Basic Principles (Kluwer Academic/Plenum Publishers, New 
York, 2003). 

C. V. Doom and D. D. Nobel, Physica, 22, 338 (1956). 



I. A. Davidson, R. T. Moug, F. Izdebski, C. Bradford, and K. A. Prior, [physica status solidi (b), 247, 1396 (2010) 



J. M. Hartmann. J. Cibert. F. Kanv , H. Mariette, M. Charleux, P. AUeysson, R. Langer, and G. Fouillct, 80, 6257 (1996) 



Smart Mater. Struct., 15, S (2006) 



D. A. Evans, A. G. McGlynn, B. M. Towlson, M. Gunn, D. Jones, T. E. Jenkins, R. Winter, and N. R. J. Poolton, Journal 
of Physics: Condensed Matter, 20, 075233 (2008). 

L. I. Berger, Semiconductor Materials (CRC Press, Boca Raton, FL, 1997). 

O. Zakharov, A. Rubio, X. Blase, M. L. Cohen, and S. G. Louie, |Phys. Rev. B, 50, 10780 (1994)| . 

While PBEsol performed well for the 18 solids studied, LSDA remained "unsurpassed" for semiconductors. See Reference**, 
page 136406-4, first paragraph. 

J . Sun, M. Ma rsman, G. I. Csonka, A. Ruzsinszky, P. Hao, Y.-S. Kim, G. Kresse, and J. P. Perdew, 



0{5117 (201i; 



Phys. Rev. B, 84 



Phys. Rev. Lett., 103, 026403 (2009)| , see also 



70 
71 
72 
73 
74 
75 
76 
77 
78 
79 



82 
83 
84 
85 
86 
87 
88 



R. Armiento and A. E. Mattsson, [Phys. Rev. B, 72, 085108 (2005) , 
J. P. Perdew, A. Ruzsinszky, G. I. Csonka, L. A. Constantin, and J. Sun, 
the Erratum*'', Phys. Rev. Lett., 106, 179902 (2011). 
P. Mori-Sanchez, A. J. Cohen, and W. Yang, Phys. Rev. Lett., 102, 066403 (2009). 

E. N. Brothers, A. F. Izmaylov, J. O. Normand, V. Barone, and G. E. Scuseria, J. Chem. Phys., 129, 011102 (2008). 
I. D. Prodan, G. E. Scuseria, and R. L. Martin, Phys. Rev. B, 76, 033103 (2007). 
L D. Prodan, G. E. Scuseria, and R. L. Martin, Phys. Rev. B, 73, 045104 (2006). 

I. D. Prodan, J. A. Sordo, K. N. Kudin, G. E. Scuseria, and R. L. Martin, J. Chem. Phys., 123, 014703 (2005). 

C. Alibert, A. JouUie, A. M. JouUie, and C. Ance, Phys. Rev. B, 27, 4946 (1983). 

H. Matthieu, D. Auvergne, P. Merle, and K. C. Rustagi, Phys. Rev. B, 12, 5846 (1975). 
A. JouUie, B. Girault, A. M. JouUie, and A. Zien-Eddine, Phys. Rev. B, 25, 7830 (1982). 
M. Grundmann, Physics of Semiconductors (Springer, Berlin, 2010). 

P. C. E. Safa Kasap, Springer Handbook of Electronic and Photonic Materials (Springer, Berlin, 2006). 

E. N. Economou, The Physics of Semiconductors Essentials and Beyond (Springer, Berlin, 2010). 
J.-W. Song, K. Yamashita, and K. Hirao, J. Chem. Phys., 135, 071103 (2011). 

M. Marsman, J. Paier, A. Stroppa, and G. Kresse, Journal of Physics: Condensed Matter, 20, 064201 (2008). 

F. El-Mellouhi, E. N. Brothers, M. J. Lucero, and G. E. Scuseria, Phys. Rev. B., 83, 205128 (2011). 
M. Snure and A. Tiwari, Appl. Phys. Lett., 91, 092123 (2007). 

D. O. Scanlon, A. Walsh, and G. W. Watson, Che m Mater., 21, 4568 (2009). 



F. Trani, J. Vidal, S. Botti, and M. A. L. Marques, Phys. Rev. B, 82, 085115 (2010) 



J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys., 124, 219906 (2006). 
J. P. Perd ew, A. Ruzsinszky, M. Ernzerhof, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and 



K. Burke, 



Phys. Rev. Lett., 102, 039902 (2009) 



J. P. Perdew, A. Ruzsinszky, G. 1. Csonka, L. A. Constantin, and J. Sun, Phys. Rev. Lett., 106, 179902 (2011). 



